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Composition, structure and stability of RuC^llO) 
as a function of oxygen pressure 
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Using density-functional theory (DFT) we calculate the Gibbs free energy to determine the lowest- 
energy structure of a RuO2(110) surface in thermodynamic equilibrium with an oxygen-rich envi- 
ronment. The traditionally assumed stoichiometric termination is only found to be favorable at low 
oxygen chemical potentials, i.e. low pressures and/or high temperatures. At realistic O pressure, 
the surface is predicted to contain additional terminal O atoms. Although this O excess defines 
a so-called polar surface, we show that the prevalent ionic model, that dismisses such termina- 
tions on electrostatic grounds, is of little validity for RuO2(110). Together with analogous results 
obtained previously at the (0001) surface of corundum-structured oxides, these findings on (110) 
rutile indicate that the stability of non-stoichiometric terminations is a more general phenomenon 
on transition metal oxide surfaces. 



I. INTRODUCTION 

Density-functional theory (DFT) is often argued to be 
a zero temperature, zero pressure technique. As such, 
the results of static total energy calculations at surfaces 
have to be transferred with considerable care to typical 
high-pressure applications like e.g. catalysis — a situa- 
tion which finds its correspondence in the experiment in 
form of the ex-situ methods of ultra high vacuum (UHV) 
surface science. Unfortunately, often enough an extrap- 
olation of the low-pressure results to technical processes 
taking place at ambient atmosphere is not possible at all, 
which has been coined with buzz words like pressure and. 
materials gap (e.g. see the discussion in Stampfl et ala 
and references therein). 

Trying to bridge these gaps, one needs to determine 
the equilibrium composition and geometry of a surface 
in contact with a given environment at finite tempera- 
ture and pressure. Then, the stable surface structure 
results from the statistical average over adsorption and 
desorption processes, and this means to employ thermo- 
dynamics. When we aim to describe experiments that 
are conducted at constant pressure and temperature, the 
appropriate thermodynamic potential to consider is the 
Gibbs free energy, G(T,p). If DFT total energies en- 
ter in a suitable way into the calculation of G(T,p) for 
a material surface, an ab initio thermodynamics results, 
which extends the predictive power of the first-principles 
technique to a more relevant temperature and pressure 
range. 

This scheme has been successfully applied to e.g. ad- 
dress rthe surface termination of corundum-structured 
oxidesaH, and we will use it here to determine the com- 
position and lowest energy structure of RuO2(110) in 
equilibrium with an oxygen atmosphere. At oxide sur- 
faces typically only so-called autocompensated, stoichio- 
metric terminations are consideredQu, i.e. they are be- 
lieved to be more favorable than the other, so-called polar 
tcrminationsQ for two reasons: First, they often involve 



a minimum of truncated bonds at the surface, and, sec- 
ond, in a purely electrostatic model, in which all oxide 
ions would be in their bulk formal oxidation state, po- 
lar surfaces would be charged and should thus exhibit an 
infinite surface energy. 

In contrast to these arguments, our ab initio thermo- 
dynamics calculations show that on RuO2(H0) a polar 
surface termination with excess oxygens is stabilized at 
higher O chemical potentials. Hence, depending on the 
experimental conditions, either the stoichiometric UHV 
or the hitherto unaccounted high-pressure termination 
can be present, which must be considered in the mod- 
elling of physical processes occuring at this surface, as e.g. 
catalytic reactions. We will also show that the rejection 
of polar surfaces on electrostatic grounds is not valid, as 
the strong dipole moment can be considerably reduced by 
surface relaxation and electron rearrangement: Rather 
than imagining the surface as created by a truncation of 
the bulk-stacking sequence at a certain bulk-like plane 
(which is the basis of the electrostatic divergence argu- 
ment), one should instead view the surface as a new ma- 
terial, in which the structural and electronic degrees of 
freedom of the top atomic layers allow a significant mod- 
ification of the properties of the atoms in the bulk. 

These results for the rutile-structured RuO2(110) are 
analogous to previous findings at the (0001) surface of 
corundum-structured oxideairtj, supporting the argument 
that polar terminations, particularly at realistic pressure, 
are a more general phenomenon on transition metal oxide 
surfaces. 



II. THEORETICAL 

This Section describes the thermodynamic formalism 
and how it is combined with DFT total energy calcula- 
tions. For the sake of clarity this is all referred explicitly 
to the present application to RuO2(110) in an oxygen 
atmosphere. However, the generalization to other com- 



pounds, M x O y , and even to an environment, that con- 
tains multiple gas phase species and not just oxygen, is 
obvious. 



A. The surface free energy 

We consider a surface in contact with an oxygen at- 
mosphere described by an oxygen pressure, p, and tem- 
perature, T. This means that the environment acts as 
a reservoir, because it can give (or take) any amount of 
oxygen to (or from) the sample without changing the 
temperature or pressure. The appropriate thermody- 
namic potential required to describe such a system is the 
Gibbs free energy, G(T,p,N Ru ,No), which depends also 
on the number of Ru, N Ru , and O, No, atoms in the sam- 
ple. The most stable surface composition and geometry 
is then the one that minimizes the surface free energy, 
j(T,p), defined as 



l{T,p) 



1 
A 



[G(T,p,N na ,N )- (1) 

N Ru fi Ru (T,p) - N Ho(T,p) ]. 



Here, /i Ru and po ar e the chemical potentials of a Ru 
and an O atom, respectively, and ~/(T,p) is normalized 
to energy per unit area by dividing through the surface 
area, A. 

If the surface system is modeled by a slab with two 
equivalent surfaces, eq. (Q) reads 



l{T,p) = — [ G slah (T,p,N Ru ,N )- 

- N R ufi Ru (T,p) - N n (T,p) 



(2) 



Now, A is the area of the surface unit cell and N Ru and 
A^o are the numbers of Ru and O atoms in the three- 
dimensional super cell. 

In eq. (0) the chemical potentials of O and Ru enter 
in a symmetric way. However, if there is enough bulk 
material, so that it acts as a thermodynamic reservoir, 
the potentials are in fact no longer independent, but are 
related by the Gibbs free energy of the bulk oxide 



MRu (T,p) + 2fi (T,p) = g$$ 2 (T,p) 



(3) 



where lower case g is henceforth used to denote a Gibbs 
free energy per formula unit. Inserting this constraint 
into eq. (0) leads to 



7(T,J0 = i 



G slab (T,p,N Rn ,N ) 



(4) 



WruS&(T,p) + (2N Ru -N )no(T,p) 



which shows how the surface free energy depends now 
only on the oxygen chemical potential. 



B. The range of allowed O chemical potentials 

As suggested by experimental conditions, we like to use 
the oxygen chemical potential to discuss the dependence 
of the surface on the O2 pressure and temperature. It 
is important to note, that experimentally (and assuming 
that thermodynamic equilibrium applies) /io cannot be 
varied without bounds. If it gets too low, all oxygen 
would leave the sample, i.e. the oxide would decompose 
into solid Ru and oxygen gas, which would start with the 
formation of Ru crystallites at the surface. Thus, 



max(» R u(T,p)) = g R T(T, P ), 



(5) 



where ff R u' k (T, p) is the Gibbs free energy of metallic 
ruthenium. Together with eq. (0) and using the T = K 
and p = atm limit for the bulk energies, we will employ 

mm( M o(T,p)) I i(0&(O,O) - 5 R u u lk (0,0) ) , (6) 

to mark the "oxygen-poor limit" (or equivalently "Ru- 
rich" limit) in the graphs discussed below. This is a good 
estimate of the real limit and most importantly, it is a 
theoretically well defined reference-point on the /io-axis. 
On the other hand, the most oxygen-rich conditions 
can be defined as the point, beyond which gas phase O 
would start to condensate on the sample. However, in 
the temperature and pressure range we are interested in 
a condensed 02-solid phase does not exist (the critical 
temperature of O2, i.e. above which gas and liquid phase 
are degenerate, is T c sw 150 K). Thus, similarly as above, 
an appropriate and well-defined estimate of the upper 
limit of the oxygen chemical potential is 



;(Mo(T,p)) = l/2££f\ 



(7) 



where _EQ tal is the total energy of a free, isolated O2 
molecule at T = K. 

Then, introducing the Gibbs free energy of formation, 
AG f (T,p), of the oxide, 

AG f (T,p) = g R ^ 2 ( T >P)-9Z k (T,p)-g^(T,p), (8) 

where g as (T,p) is the Gibbs free energy of an O2 
molecule, we see that the range of oxygen chemical po- 
tentials between our theoretical tick marks is 

iAG/(0,0) < ^ {T,p)-\e^ 1 < 0. (9) 

We compute AG/(0,0) = —3.16 eV per formula unit, 
which compares very well with the experimental Gibbs 
free energy of formation at standard pressure in the limit 
of low temperatures, AG° f (T -» OK, latm) = -3.19 cV 
per formula unitQ. 

It is important to notice, that our "tick marks" for the 
oxygen-rich and oxygen-poor conditions are theoretically 
well defined limits, yet they only represent an estimate of 



the truly accessible range of the oxygen chemical poten- 
tial. The range between our tick marks is 1/2 AG/(0, 0), 
but in reality the accessible range is 1/2 AG/ (T, p) , i.e. it 
is temperature and pressure dependent. At T — 1000 K 
and p = 1 atm the Gibbs free energy of formation has 
increased by 0.63 eV compared to the aforementioned 
T — > K valueQ. Keeping the ensuing variation of 0.3 eV 
in the accessible range of potentials in mind, we will al- 
ways show the resulting curves in our below discussed 
figures also some tenths of an eV outside the "oxygen- 
rich" and "oxygen-poor" tick marks. 



C. The oxygen-poor limit as a safe reference 

Total energies for extended systems are typically more 
accurately described by actual DFT calculations than 
those for atoms and molecules. With respect to eq. (0), 
it is therefore suitable to rewrite it to 



~ (ffR U u0 2 (0, 0) - <&f (0, 0)) < mo (T, p) < 

^(sruO 2 (0,0) -0&f (0,0)) + ^AG/(0,0). 

If we then insert the oxygen-poor limit into eq. 
obtain for the surface free energy 



(10) 



— poor(-i iP) — 


1 

2A 


G slah (T,p,N Ru ,N ) - N Ru gZo 2 ( T >P) 




(Vru-^) « lk (0,0) - 5 S 2 (0,0)) 



(11) 



Likewise, the oxygen-rich limit turns out to be 



7o- 



ch(T,p) =7o^poor(T,p) 



(12) 



^(tfR»-^)AG,(.UM. 



The result of this rewriting of eq. (0) is that atomic 
or molecular quantities do not enter into the calculation 
of the oxygen-poor limit, i.e. eq. (O), at all, which 
thus defines a safe reference involving only bulk or slab 
quantities. 

On the other hand, AG/ (0,0) depends on the O2 to- 
tal energy, and AG/ (0,0) defines the slope of the lines 
representing the surface free energy as a function of /10 : 
The slope is 



1 
2A 



N Ru 



Nr 



AG/ (0,0), 



(13) 



cf. eq. (|1 2[), and sometimes AG/ (0,0) may be affected 
by the error in E ° tal , in which case it might be preferable 
to use its experimental value. Yet, for the present case of 
Ru0 2 , we note that our DFT-GGA result for the Gibbs 
free energy of formation is very close to the experimental 



value (s. above). Thus, here the error in E 0tal , which 
clearly exists, cancels out and the calculated slopes are 
very accurate. As a consequence, and in contrast to com- 
mon belief, we note that the bulk total energy of Ru02 

must therefore have a similar error as £f? tal - otherwise 

2 
the apparent error cancelation in AG/ (0,0) would not 

occur. 

We finally note in passing, that eqs. ( |ll| ) and ( |l2| ) 
nicely reflect the physics behind the dependence on the 
O chemical potential: While a stoichiometrically termi- 
nated surface structure (N Ru = No/2) will exhibit a 
constant surface free energy as a function of no(T,p), 
a termination with an O excess (deficiency) will become 
more and more favorable (unfavorable) with increasing 
Ho(T,p), i.e. higher O pressure and/or lower tempera- 
ture. 



D. Gibbs free energies vs. total energies 

The formalism as described up to this point is entirely 
based on the Gibbs free energies of the system, whereas 
we intend to provide as input total energies from DFT 
calculations. Therefore, we will now outline how both 
quantities are related, and under which approximations 
they might be equated to each other. 

DFT total energies are evaluated for a certain volume, 
V, of the unit cell. The resulting E total (V,N Rn ,N ) 
is related to a thermodynamical quantity only in a re- 
stricted way: It corresponds to the Helmholtz free energy 
at zero temperature and neglecting zero-point vibrations. 
In general, the Helmholtz free energy can thus be written 
as 



F(T, V, A Ru , N ) = E tota \V, N Ru , N Q ) 



F vih -(T,V,N Ru ,N ), 



with 



TS vih -(T,V,N Rn ,N ) 



(14) 



(15) 



comprising all contributions, which depend on vibra- 
tional modes in the system. Here, E vlh - and S vlb - are 
the vibrational energy (including the zero-point energy) 
and entropy respectively. In turn, the Helmholtz free 
energy is associated to the Gibbs free energy via 



G{T,p,N Ru , N ) 



F(T,p,N Ru ,N ) - 
pV(T :P ,N Rn ,N ). 



(16) 



Checking first on the pV term, we find from a simple 
dimensional analysis, that its contribution to the sur- 
face free energy (normalized to the surface area) will be 
[pV/A] = atm A 3 / A 2 ~ 10~ 8 meV/A 2 . As we are only 
interested in a pressure range that will not exceed about 
100 atm, this contribution is completely negligible com- 
pared to the Helmholtz free energy, which is of the order 
of me V/ A 2 . 



This leaves as the only additional contribution to 
G(T, p, Nr u ,No) apart from the DFT total energy only 
the vibrational term F vih (T, V, N Ru , No). Using the 
phonon density of state (DOS), cr(u>), this vibrational 
part of the free energy can be written as an integral over 
the modes, u>, 

F vih -(T,V,N Ru ,N ) = [duF vih iT,Lj)a(u;), (17) 



10.0 



where an analytic expression for _F vlb - (T, u>) is given in 
the appendix. 

Plugging this into eq. (pi) we obtain for the vibra- 
tional contribution to the surface free energy of a sto- 
ichiometric termination (Ar u = No/2) at the O-poor 
limit. 



'O — poor 



1 

2A 



(T,V) = 

da; F v[h -(T,Lu) 



(18) 



slab , 



u;) - N Ru aZo 2 H) 



To get an estimate of its value, we use the Einstein model 
and approximate the phonon DOS by just one character- 
istic frequency for each atom type. If we further consider, 
that the vibrational mode of the topmost layer of Ru and 
of O might be significantly changed at the surface, we 
thus have d>Q Ulk and w{^" lk as characteristic frequencies of 
O and Ru in Ru02 bulk, as well as tj™ rf - and wf^ as 
the respective modes at the surface. With this simplified 
phonon DOS, eq. ( |l8| ) reduces to 



J V 6- P oor(T,V) 

3 
2A 



(19) 



TpVlb.frri -SUrf.N T^Vlb.frri -bulloN ■ 

F vib -(T, cDg 11 *-) - F vih {T,u! nlk )) 



Hence, we see that in the O-poor limit 7o— poorC^P) °^ 

a stoichiometric termination arises primarily out of the 
difference of the vibrational modes at the oxide surface 
with respect to their bulk value. To arrive,a.t numbers, 



we use u), 



bulk 



80 meV and iD R ^ k 



25 meVB'u, and allow 
a 50% variation of these values at the surface, to plot eq. 
( |l9| ) in Fig. [y in the temperature range of interest to our 
study. As these particular values for the characteristic 
frequencies are not well justified, but should only be con- 
sidered as rough estimates, we also include in the graph 
the corresponding To^- r (T,p), ^ these values changed 
by ±50%. From Fig. [Owe see that the vibrational contri- 
bution to the surface free energy stays within 10 meV/A 2 
in all of the considered cases, and that the uncertainty 
in the characteristic frequencies translates primarily to 
variations of Jo-vootC^'P) a t l° w temperatures, where 
the value of the latter is very small anyway. 

We have also computed the vibrational contribution 
to the surface free energy of non-stoichiomctric termina- 
tions in an analogous manner. There, the formulas turn 




1000 



Temperature (K) 

FIG. 1. Vibrational contribution to the surface free energy, 
of a stoichiometric termination, cf. eq. (119), in the temper- 
ature range of interest to the present study. The Ru and O 
modes are approximated in the Einstein model by character- 
istic frequencies, cD ulk = 80 meV and w R u lk =25 meV. Shown 
is the contribution, if the vibrational modes at the surface 
differed by ±50% from these bulk values (solid lines). To 
assess the dependence on the value chosen for the character- 
istic frequencies, the latter are varied by ±50% (dashed and 
dotted lines respectively). In all cases, the vibrational con- 
tribution stays below 10 meV/A 2 in the whole temperature 
range considered. 



out more lengthy, cf. eq. (Jill) , and the vibrational con- 
tribution includes not only differences between bulk and 
surface vibrational modes, but also absolute F vlb (T, u) 
terms due to the excess or deficient atoms. Yet, even then 
the vibrational contribution stays within ±10 meV/A 2 
similar to the above described stoichiometric case. In 
conclusion, we therefore take this value to represent a 
good upper bound for the vibrational influence on the 
surface free energy. 

Such a ±10 meV/A 2 contribution is certainly not a 
completely negligible factor, yet as we will show below 
it is of the same order as the numerical uncertainty in 
our calculations. Furthermore, as will become apparent 
in the discussion of the results, this uncertainty does not 
affect any of the physical conclusions drawn in the present 
application. Hence, we will henceforth neglect the com- 
plete vibrational contribution to the Gibbs free energy, 
leaving only the total energies E total (V, A Ru , No) as the 
predominant term. In turn, this allows us to rewrite eq. 
© to 



70-poor(r,p) 



1 

2A 



£ slab (T/,A Ru ,A ) 



(20) 



, « 2 W-(^Ru-^ 



KT(V) 



which now contains exclusively terms directly accessible 
to a DFT calculation. We stress, that this approximation 
is well justified in the present case, but it is not a gen- 



TABLE I. no(T,p°) in the temperature range of inter- 
est to our study. The entropy and enthalpy changes used 
to obtain fio(T,p ) via eq. (E2J) ares-taken from the JANAF 
thermochemical tables at p° — latmlij. 



F. DFT basis set and convergence 



T 


Mo(T,p°) 


T 


Mo(T,P°) 


100 K 


-0.08 eV 


600 K 


-0.61 eV 


200 K 


-0.17 eV 


700 K 


-0.73 eV 


300 K 


-0.27 eV 


800 K 


-0.85 eV 


400 K 


-0.38 eV 


900 K 


-0.98 eV 


500 K 


-0.50 eV 


1000 K 


-1.10 eV 



eral result: There might well be applications, where the 
inclusion of vibrational effects on the surface free energy 
can be crucial. 



E. Pressure and temperature dependence of /xo (T, p) 

Having completely described the recipe of how to ob- 
tain j(T,p) as a function of the O chemical potential, the 
remaining task is to relate the latter to a given temper- 
ature, T, and pressure, p. As the surrounding O2 atmo- 
sphere forms an ideal gas like reservoir, we can obtain in 
the appendix the following expression 



fJo{T,p) = p, {T,p°) 



1/2 kTlnl — 
\p° 



(21) 



which already gives us the temperature and pressure de- 
pendence, if we only know the temperature dependence 
of no(T,p°) at one particular pressure, p° . 

We choose as the zero reference state of ^o(T,p) 
the total energy of oxygen in an isolated molecule, i.e. 
Ho(0K,p) = 1/2 £&° tal = 0. With respect to this zero, 
reached at our theoretical oxygen- rich limit, cf. eq. (|9j), 
Ho(T,p°) is then given by 



Vo(T,p°) 



M g- rich (0K,p°) 



1/2AG(A7>°,0 2 ) = 
1/2 (H(T,p°,O 2 )-H(0K,p°,O 2 )) - 
1/2 T (S(T,p°,0 2 ) - S(0K, P °, 2 )) , (22) 



where we have used the relation G = H — TS between 
the Gibbs free energy and the enthalphy, H . This allows 
us to obtain the aspired temperature dependence simply 
from the differences in the enthalpy and entropy of an 2 
molecule with respect to the T = K limit. For standard 
pressure, p° = 1 atm, these values axe e.g. tabulated 
in the JANAF thermochemical tablesEHl. Plugging them 
into eq. ( B2 ) leads finally to (j,q(T,p°), which we list in 
Table I 

Together with eq. (|l]) the O chemical potential can 
thus be obtained for any given (T,p) pair. Although we 
prefer to conveniently present the resulting surface en- 
ergies as a one-dimensional function of ^0 {T, p) , we will 
often convert this dependence also into a temperature 
(pressure) dependence at a fixed pressure (temperature) 
in a second x-axis to elucidate the physical meaning be- 
hind the obtained curves. 



The DFT input to eq. (20) has been obtained using the 
full-potenfi,a| linear augmented plane wave method (FP- 



LAPW)tirt3 within the generalized gradient approxima- 
tion (GGA) of the exchange-correlation functionafJ. For 
the RuO2(H0) surface calculation we use a symmetric 
slab consisting of 3 rutile 0-(RuO)-0 trilayers, where all 
atomic positions within the outermost trilayer were fully 
relaxed. A vacuum region of « 11 A is employed to de- 
couple the surfaces of consecutive slabs in the supercell 
approach. Test calculations with 5 and 7 trilayered slabs, 
as well as with a vacuum region up to « 28 A confirmed 
the good convergence of this chosen setup with variations 
of j(T,p) smaller than ± 3 meV/A 2 . Allowing a relax- 
ation of deeper surface layers in the thicker slabs did not 
result in a significant variation of the respective atomic 
positions, nor did it influence the near-surface geome- 
try as obtained in the calculations with the standard 3 
trilayer slabs. To ensure maximum consistency, the cor- 
responding Ru0 2 bulk computations are done in exactly 
the same (110) oriented unit cell as used for the slabs, 
in which the prior vacuum region is simply replaced by 
additional Ru0 2 trilayers. 

The FP-LAPW basis set is taken as follows: R^ T =1.8 
bohr, i?MT =1-3 bohr, wave function expansion inside 
the muffin tins up to /J^ x = 12, potential expansion up 
to l^ x = 4. For the RuO 2 (110) slabs the Brillouin zone 
integration was performed using a (5 x 10 x 1) Monkhorst- 
Pack grid with 15 k-points in the irreducible part. The 
energy cutoff for the plane wave representation in the in- 
terstitial region between the muffin tin spheres was 17 
Ry for the wave functions and 169 Ry for the potential. 
Checking on the convergence, the surface free energies 
of the three possible (1 x 1) RuO 2 (110) truncations dis- 
cussed below were found unchanged to within 1 meV/A 2 
by increasing the k-mesh to a (7 x 14 x 1) Monkhorst- 
Pack grid with 28 k-points in the irreducible part. A 
larger interstitial cutoff up to 24 Ry reduced the abso- 
lute values of the three j(T,p) by up to 10 meV/A 2 , yet 
as all of them were reduced alike, their respective differ- 
ences (which are the only relevant quantities entering the 
physical argument) stayed constant to within 5 meV/A 2 . 

In total we thus find the numerical accuracy of the 
calculated surface free energies with respect to the su- 
percell approach and the finite basis set to be within 10 
meV/A 2 , which will not affect any of the physical conclu- 
sions drawn. Note, that the stated imprecision does not 
include possible errors introduced by more general defi- 
ciencies of the approach, namely the use of the GGA as 
exchange-correlation functional, on which we will com- 
ment below. 



III. RESULTS 



a) RuO 2 (110) - O bridge 
Qbridge Ru cus,5f Ru bridge,6f 




b) RuO 2 (110) - O 

n CUS,6f qCUS 



cus 




c) RuO 2 (110) - Ru 

Ru cus,5f Ru bridge,4f 




FIG. 2. Three possible terminating planes of the rutile 
(110) surface: a) Stoichiometric RuO 2 (110)-O bridgc termi- 
nation with five-, six- and twofold coordinated Ru cus,5f , 

Ru bridgc,6f and Qbridge atomg respectively . b ) R U 2 (H0)-O CUS 

termination, where additional O cus atoms sit atop the for- 



merly undercoordinated Ru c 



termination, which lacks the O rl ge 
the stoichiometric termination (Ru - 
— small, dark spheres). 



atoms, c) RuO 2 (110)-Ru 



atoms in comparison to 
= large, light spheres, O 



A. RuO2(110) surface structure 

RuC>2 crystallizes in the rutile bulk structure, in which 
every metal atom is coordinated to six oxygens, and every 
oxygen to three metal neighborstJ. The oxygens forming 
an octahedron around each Ru atom are not all equiv- 
alent, but can be distinguished into four basal and two 
apical O atoms with calculated O-Ru bondlengths of 2.00 
A and 1.96 A respectively. We notice that along the 
(110) direction this structure can then be viewed as a 
stacking sequence of 0-(RuO)-0 trilayers, in which each 
trilayer is simply composed of an alternating sequence 



of in-plane and up-right oriented oxygen-ruthenium co- 
ordination octahedra, cf. Fig. gb. Cut along (110), the 
rutile bulk structure can therefore exhibit three distinct 
terminations of (1 x 1) periodicity, depending at which 
plane the trilayer is truncated, cf. Fig. Ha-c. 

Traditionally, the stoichiometric RuO 2 (110)-O bridgc 
termination is believed to be the most stable one for all 
(110) surfaces of rutile-structured crystalsQiJ, as in the 
ionic model it leads to an uncharged surface and cuts 
the least number of bonds: While the Ru brldgc ' 6f atoms 
possess their ideal sixfold O coordination with two of 
their basal oxygens forming the terminal O brldgc atoms, 
only the Ru cus,5f lack one apical on-top O, cf. Fig. ||a. 
Note, that we will use a nomenclature for the surface Ru 
atoms, where apart from a site specific characterization 
(e.g. cus for the coordinatively unsaturated site in the 
stoichiometric termination) also the number of direct O 
neighbors (e.g. 5f for 5-fold coordination) is stated. Vice 
versa, we indicate for the surface O atoms to which spe- 
cific site they bind (e.g. O bridgo binds to the Ru bridgc > 6f 
atoms). 

Alternatively, in the second possible RuO2(110)-O cus 
termination, cf. Fig. gb, terminal O cus atoms oc- 
cupy sites atop of the formerly undercoordinated Ru cus ' 
atoms, so that now all metal atoms in the surface pos- 
sess their ideal sixfold coordination. This is, of course, 
paid for by the presence of both the only twofold and 
onefold coordinated O brldge and O cus atoms, respectively. 
Finally, the third RuO2(H0)-Ru termination exhibiting 
the mixed (RuO) plane at the surface is achieved by re- 
moving the O brldgo atoms from the stoichiometric termi- 
nation, cf. Fig. |p. Here, no undercoordinated oxygens 
are present any more, yet, at the expense of the four- and 
fivefold bonded R u b "dgc,4f and Ru cus ' 5f atoms. 



B. Prediction of a high pressure termination 

The calculated surface free energies of the three pos- 
sible terminations are shown in Fig. pi As ex- 
plained in connection with eq. (|l2|), the RuOa(HO)- 
qcus (R u O 2 (110)-Ru) termination with an excess (defi- 
ciency) of O at the surface becomes more and more fa- 
vorable (unfavorable) towards the O-rich limit, while the 
stoichiometric RuO2(110)-O brldgc termination exhibits a 
constant j(T,p). Indeed, we find the traditionally as- 
sumed stoichiometric RuO2(H0)-O brldgc surface to be 
the most stable over quite a range of oxygen chemical 
potentials above the O-poor limit. To make this range a 
bit more lucid, we have used eq. (^lj) to convert no(T,p) 
into a physically more comprehensible pressure depen- 
dence at a fixed temperature, cf. the top x-axis of Fig. 

The temperature of T — 600 K corresponds to a typ- 
ical annealiag-,tenjiperature employed experimentally for 
this systerro'Ej~t£l. From the resulting pressure scale we 
see, that the stability of the stoichiometric termination 
extends therefore roughly around the pressure range cor- 
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FIG. 3. Surface free energies, y(T,p), of the three 
RuC>2(110) terminations depicted in Fig. bl Addition- 
ally shown as dashed line is the surface free energy of a 
RuO2(110)-O cus termination, in which only every second O cus 
site along the trenches is occupied. The dotted vertical lines 
indicate the allowed range of the oxygen chemical potential, 
Ho{T,p), using 1/2 Eq^ 1 as zero reference as explained in 
Section IIB. In the top a;-axis, the dependence on ^o(T,p) 
has been cast into a pressure scale at a fixed temperature of 
T = 600 K. 



responding to UHV conditions. 

Yet, this is different at higher O pressures, where the 
RuO2(110)-O cus termination becomes the most stable 
surface structure, cf. Fig. ||. In the O-rich limit, it 
exhibits a 7°~ rlch (T,p), which is by 49 meV/A 2 lower 
than the one of the stoichiometric RuO2(110)-O brldge sur- 
face, i.e. the deduced crossover between the two termi- 
nations is far beyond the estimated uncertainty of ~ ±10 
meV/A 2 due to the neglection of the vibrational contri- 
bution to the Gibbs free energies and due to the finite 
basis set. This estimate docs, however, not comprise 
the more general error due to the use of the GGA as 
exchange-correlation functional. To this end, we have 
also calculated the surface free energies of the two com- 
peting terminations within the local density approxima- 
tion (LDA)E3. Although the absolute values of both 
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(T,p) turn out to be by « 15 meV/A 2 higher, 
their respective difference is almost unchanged, which is 
eventually what determines the crossover point of the two 
lines in Fig. |. The RuO 2 (H0)-O cus termination is there- 
fore the lowest energy structure for no(T,p) > —0.85 eV 
in the LDA, which is almost the same as the [io(T,p) > 
—0.93 eV found with the GGA, shown in Fig. a. Con- 
sequently, while the choice of the exchange-correlation 
functional may affect the exact transition temperature 
or pressure, the transition per se is untouched. In turn, 
we may safely predict the stability of a polar surface ter- 
mination on RuO2(110) at high O chemical potential, 
corresponding e.g. to the pressure range typical for cat- 
alytic applications. 



Note, that Fig. || summarizes only the j(T,p) of the 
three (lxl) terminations, which arise by truncating the 
Ru02 crystal at bulk-like planes in the (110) orientation. 
Yet, it is a priori not clear, that the terminal atoms at 
the surface must be in the sites corresponding to the bulk 
stacking sequence. To check this, we have additionally 
calculated the surface free energies of surface structures, 
where the O brldge (O cus ) atoms occupy atop (bridge) sites 
over the R A1 brldgG ' 5f (Ru cus,6f ) atoms, instead of their nor- 
mal bridging (atop) configuration. In both cases, we find 
the j(T,p) considerably higher, which excludes the pos- 
sibility that these adatoms occupy non bulk-like sites at 
the surface. Similarly, a RuO2(110)-O cus termination, 
where the O brldgc atoms have been removed, can also 
safely be ruled out as alternative for a stoichiometrically 
terminated surface. 



C. Lateral interaction between O cus atoms and 
vacancy concentration 

The reasoning of the last Section leaves in fact only 
the RuO 2 (H0)-O bridgc and RuO 2 (H0)-O cus termination 
as the relevant surface structures stabilized in UHV and 
under high O pressure respectively. Both differ from 
each other only by the presence of the additional O cus 
atoms, which continue the bulk stacking sequence by fill- 
ing the vacant sites atop of the formerly undercoordi- 
nated Ru cus ' 5f atoms. The way the RuO 2 (H0)-O cus ter- 
mination is formed at an increasing oxygen chemical po- 
tential will therefore depend significantly on the details 
of the lateral interaction among this adatom species. 

The RuO 2 (110)-O bridgo surface has a trench-like struc- 
ture with a distance of 6.43 A between the rows formed 
by the O brldgc atoms, cf. Fig. ||a. This renders any lat- 
eral interaction between O cus atoms adsorbed in neigh- 
boring trenches rather unlikely. On the other hand, the 
distance between two O cus atoms occupying neighbor- 
ing sites along one trench is only 3.12 A. To check on 
the corresponding interaction we calculated the surface 
free energy of a RuO2(110)-O cus termination in a (2 x 1) 
supercell, in which the O cus atoms occupied only every 
other site along the trenches. The corresponding 7(T,p) 
is drawn as a dashed line in Fig. ||. As now only half of 
the excess O cus atoms are present, the slope of this curve 
has to be one half of the slope of the line representing 
the normal RuO2(H0)-O cus termination, cf. eq. (|l2|). 

Interestingly, both curves cross the stoichiometric 
RuO 2 (H0)-O bridgc line in exactly the same point. This 
can only be understood by assuming a negligible lateral 
interaction between neighboring O cus atoms: If there was 
an attractive (repulsive) interaction between them, then 
it would be favorable (unfavorable) to put O cus atoms 
as close to each other as possible. In turn, the (2x1) 
overlayer of O cus atoms, in which only every other site 
is occupied, would be less (more) stable than the nor- 
mal RuO2(110)-O cus termination, where all neighboring 



sites are full. Consequently, the stability with respect to 
the stoichiometric termination would be decreased (en- 
hanced), leading to a later (earlier) crossover point in 
Fig. [| That both calculated lines cross the RuO2(110)- 
Qbndgo jj ne a £ -j-^g same p i n t i s therefore a reflection of 

a negligible lateral interaction between the O cus atoms. 

Additionally, we compute a very high barrier of almost 
1.5 eV for diffusion of O cus atoms along the trenches, indi- 
cating that the latter species will be practically immobile 
in the temperature range, where the oxide is stable. This 
together with the small lateral interaction indicates that 
at increasing O chemical potential the RuO2(110)-O cus 
surface is formed from the stoichiometric termination by 
a random occupation of O cus sites, until eventually the 
whole surface is covered. Not withstanding, at finite tem- 
peratures there will still be a certain vacancy concentra- 
tion even at O chemical potentials above the crossover 
point of the two terminations. As the undercoordinated 
Ru cus ' atoms exposed at such a vacant O cus site, cf. Fig. 
0, might be chemically active sites for surface reactionsli3, 
it is interesting to estimate how many of these sites will 
be present under given (T,p) conditions. 

Since we have shown that each O cus site at the surface 
is filled independently from the others, we can calculate 
its occupation probability within a simple two level sys- 
tem (site occupied or vacant) in contact with a heat bath. 
The vacancy concentration follows then from a canonic 
distribution, where the energy of the two levels is given by 
the j(T,p) of the two terminations at the chosen chemical 
potential. As an example, we first address room temper- 
ature, where the RuO2(110)-O cus termination becomes 
stable at pressures higher than p ~ 10~ 22 atm. A va- 
cancy concentration of only 1% is in turn already reached 
at p ~ 10~ 17 atm, so that at this temperature there will 
only be a negligible number of vacancies on the O cus - 
terminated surface for any realistic pressure. 

However, this situation becomes completely different 
at elevated temperatures. At T = 800 K, the crossover 
to the RuO2(110)-O cus -termination occurs at p ~ 10 _1 
atm with a 10% vacancy concentration still present at 
10 2 atm. In the range of atmospheric pressures the 
RuO2(110) surface will therefore exhibit a considerable 
number of vacancies, which could explain, the , high cat- 
alytic activity reported for this materiafliSEilO ~E-X How- 
ever, although we have deliberately chosen T = 800 K as 
a typical catalytic temperature, where e.g. a maximum 
conversion rate for the_CO/C02 oxidation reaction over 
RuO2(110) was foundE3, we immediately stress that our 
reasoning is at the moment only based on the O pressure 
alone and therefore not directly applicable to catalysis 
experiments, which may also depend on the partial pres- 
sure of other reactants in the gas phase. 
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FIG. 4. Surface free energies, j(T,p), of the stoichiometric 
RuO 2 (110)-O bridgc and polar RuO 2 (110)-O cus termination. 
Shown is the effect of relaxation at the surface, with solid 
lines indicating fully relaxed surface structures and dashed 
lines the corresponding bulk-truncated geometries. The de- 
pendence on the oxygen chemical potential has been trans- 
lated into a temperature scale at 10~ 12 atm (bottom :r-axis) 
and 1 atm (top :r-axis) pressure. 



D. On the stability of polar surfaces 

As already mentioned in the introduction, the pre- 
dicted high pressure RuO2(H0)-O cus termination is tra- 
ditionally not expected as it forms a so-called polar 
surface, -which should not be stable on electrostatic 
groundsu'U. The corresponding argument is based on the 
ionic model of oxides, in which every atom in the solid is 
assumed to be in its bulk formal oxidation state. Along a 
particular direction, z, the crystal may then be viewed as 
a stack of planes with charge q, each of which contribute 
with V(z) oc qz to the total electrostatic potential. As 
this contribution diverges at infinite distances the crys- 
tal as a whole can in turn only be stable, if constructed 
as a neutral block, in which all infinities due to the in- 
dividual planes cancel. For RuO2(110), which is a type 
2 surface in Tasker's widely used classification schemeEa, 
the only neutral repeat unit is a symmetric O-(RuO)- 
O trilayer with a (-2)-(+4)-(-2) charge sequence, so that 
the only surface termination without net dipole moment 
would correspondingly be the stoichiometric RuO2(H0)- 
bridgc termination. 

On the other hand, the RuO2(110)-O cus termination 
with its extra unmatched (-2) charge plane formed by 
the O cus atoms would lead to a diverging potential and 
should thus not be stable. That we indeed find this sur- 
face stabilized at higher O chemical potentials points to 
the most obvious shortcoming of the electrostatic model, 
namely the assumption that all atoms in the solid are 
identical, i.e. that also all surface atoms are both struc- 
turally as well as electronically in a bulk-like state. In 
how much already the additional structural degrees of 



freedom at the surface influence the stability is exem- 
plified in Fig. [|, where the surface free energies of the 
two relevant terminations are compared in either a bulk- 
truncated or the fully relaxed geometry. While the small 
relaxation of the RuO2(110)-O brldgc termination hardly 
affects the y(T,p), the bulk-truncated RuO 2 (110)-O cus 
surface turns out considerably less stable compared to its 
relaxed counterpart. The by +2.5 eV strongly increased 
workfunction of the bulk-truncated RuO2(110)-O cus sur- 
face with respect to the stoichiometric termination shows 
that the addition of the O cus atoms indeed induces a con- 
siderable dipolc moment as suggested by the ionic model. 
Yet, just the relaxation lowers this workfunction again 
by 1 eV, indicating that the dipole moment can already 
be considerably reduced via a strongly shortened O cus - 
Ru cus ' 6f bondlength of 1.70 A (compared to the bulk-like 
1.96 A), therewith considerably stabilizing the surface. 

Yet, not only structurally, but also electronically the 
topmost layers in the RuO2(110)-O cus surface differ ap- 
preciably from their respective bulk counterparts. This 
is illustrated in Fig. || where we show the (x, y)- 
averaged potential along the (110) direction perpendicu- 
lar to the surface, z. In the bulk-truncated, stoichiomet- 
ric RuO2(110)-O brldgc termination in Fig. @a the elec- 
trostatic potential at the topmost 0-(RuO)-0 trilayer 
is still almost identical to the corresponding one in the 
deeper trilayers, thus enabling in this case a description 
of this surface in terms of bulk-like planes as assumed 
in the ionic model. On the contrary, we find a signifi- 
cant deviation of the potential in the outermost layers of 
the RuO2(110)-O cus termination, cf. Fig. pp, even for a 
bulk-truncated geometry. This difference is further en- 
hanced by the structural relaxation, cf. Fig. pb, so that 
the topmost (RuO)-O-O planes of this termination are 
certainly not well characterized by bulk properties, thus 
invalidating the electrostatic argument raised against this 
polar surface. 

Instead, we argue that the surface fringe composed by 
the topmost layers should be viewed as a new material, 
which properties might differ considerably from the bulk 
stacking sequence due to the additional structural and 
electronic degrees of freedom present at the surface. A 
similat,conclusion has previously been reached by Wang 
et aZ£m, who discussed the stability of oxygen-terminated 
polar (0001) surfaces of corundum-structured a-Fe2C>3 
and a-A^Oa. This indicates, that the traditionally dis- 
missed polar terminationsQ might indeed be a more gen- 
eral phenomenon, which existence could be a crucial in- 
gredient to understand the function of oxide surfaces un- 
der realistic environmental conditions. As particularly 
polar terminations with excess oxygen can be stabilized 
at increased O2 partial pressure in the gas phase, the dif- 
ferent properties of the latter should be taken into consid- 
eration, when modelling high pressure applications like 
catalysis. 
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FIG. 5. (x, y)-averaged Kohn-Sham effective, V c (2) 
(dashed line), and electrostatic potential, V CB (z) (solid line), 
along the (110) direction perpendicular to the surface, z. Ad- 
ditionally shown is the workfunction, $, as dotted line, a) 
Bulk-truncated RuO 2 (110)-O bridgc termination ($ = 5.7 eV), 
b) Bulk-truncated RuO 2 (110)-O cus termination ($ = 8.2 eV), 
and c) Fully relaxed RuO 2 (110)-O cus termination ($ = 7.2 
eV). The top s-axis marks the position of 0-(RuO)-0 layers 
in the crystal. 



E. Importance of experimental preparation 
conditions 

This influence of the O2 partial pressure on the sur- 
face morphology and function has recently become ap- 
parent in a number of studies addressing the reported 



high CO oxidation rates over Ru catalysts. While it was 
for a long time-believed that the active species is the 
Ru metal itselfc3c3, the decisive role played by oxide 
patches forrned-. undcr . c atalytic conditions was only re- 
cently realizedliZro£iE3. This was primarily due to the 
problem of preparing a fully oxidized surface in a con- 
trolled manner under UHV conditions. Yet, by means of 
more oxidizing carrier gases or higher O partial pressures 
it is now possible to circumvent this materials gapzB, 
enabling the detailed characterization of RuO2(110) do- 
mains formed on the modelr-iRjitoOOii-Burface with the 
techniques of surface scienceBllZnl3E2rE3. 

However, the results of the present study show that 
the surface termination of these domains changes with 
the O2 pressure, too. Not aware of this dependence 
a reaction mechanism for the CO oxidation was pre- 
viously proposed, situated solely on the stoichiometric 
RuO2(H0)-O bridgo termination that was characterized in 
the respective LEED study under UHV conditionsE£tE3. 
While it is presently not clear, to what extent the O cus 
atoms additionally pres ent u nder catalytic pressures are 
involved in the reactiorHEJ, this still highlights the del- 
icacy with which the results of UHV spectroscopies and 
post-exposure experiments have to be applied to model 
catalysis and steady-state conditions. 

Only very recently UHV studies achieved to stabilize 
the high-pressure RuO2(H0)-O cus terminajtian_iptention- 
ally by postdosing O2 at low temperaturesaOE^O. Tem- 
perature desorption spectroscopy (TDS) experiments 
found the corresponding excess 0£^, s atoms to be sta- 
ble up to about 300-550 K in UHVtl This agrees nicely 
with the calculated transition temperature of 450 ± 50 K 
at the crossover point between the two terminations for 
a pressure of 10 -12 aipj, presumably present during a 
typical TDS experimentEil, cf. the bottom x-axis of Fig. 
[|. Yet, the actual desorption temperature is of course sig- 
nificantly higher for the orders of magnitude higher O2 
partial pressures present in catalytic applications. This 
is exemplified by the temperature scale on the top x-axis 
of Fig. I at a pressure of p = 1 atm, representative 
for the early high-pressure experiments addressing the 
high CO/CO2 conversion rates of Ru catalystsEHa. The 
corresponding elevated transition temperature of 900 K 
shows that O cus atoms were most probably present on 
oxidized RuO2(110) domains in all of these experiments. 

The presence of the hitherto unaccounted high- 
pressure termination of RuO2(110) might therefore be 
an essential keystone to understand the data obtained 
from grown RuO2(H0) films or oxidized Ru(0001) sur- 
faces, which are unanimously prepared under highly O- 
rich conditions —We already exemplified this in a preced- 
ing publications by suggesting that the long discussed, 
largely shifted satellite peak in Ru 3d x-ray pholaemis- 
sion spectroscopy (XPS) data from such surfacesE£lt3, is 
due to the Ru cus - 6f atoms in the RuO 2 (110)-O cus termi- 
nation, which experience a significantly different environ- 
ment due to the aforementioned, very short bondlcngth 
to the O cus atoms. 



If the high-pressure termination created during the 
preparation of the crystal is frozen in during the trans- 
fer to UHV depends then on the details of the trans- 
fer itself, e.g. on such nitty-gritties of whether or not 
the temperature is kept constant while the pressure goes 
down to its base value after exposure. A dependence 
of the TDS data of an oxygen-rich Ru(0001) surface on 
these parameters has aheadv-been reported and consid- 
ered by Bottcher and NiehualJ, while the final annealing 
step to 600 K after transfer to UHV employed in the 
LEED wtpk-identifying RuO2(110) domains on oxidized 
Ru(0001)t§tL3 explains, why there only the stoichiometric 
RuO2(H0)-O bndgc termination could be characterized, 
cf. Fig. |. 

Such dependences on the experimental preparation 
have hitherto often been neglected, entailing only a low 
comparability of data sets obtained in different groups. 
Instead, the present results demonstrate that systematic 
investigations in the whole (T,p) range are required to 
fully identify the surface structure and composition of 
oxide surfaces at realistic conditions, which in turn is 
a prerequisite before tackling the long-term goal of un- 
derstanding the function of the latter in the wealth of 
everyday applications. 



IV. SUMMARY 

In conclusion we have combined density-functional the- 
ory (DFT) and classical thermodynamics to determine 
the lowest energy structure of an oxide surface in equi- 
librium with an O environment. The formalism is applied 
to RuO2(H0) showing that apart from the expected sto- 
ichiometric surface, a so-called polar termination with 
an excess of oxygen (O cus ) is stabilized at high O chem- 
ical potentials. Depending on the details of the experi- 
mental preparation conditions, either of the two termina- 
tions can therefore be present, which different properties 
have to be taken into account, when trying to understand 
the obtained data or aiming to extrapolate the results of 
UHV ex situ techniques to high-pressure applications like 
oxidation catalysis. 

A polar termination is traditionally not conceived to 
be stable within the framework of electrostatic arguments 
based on the ionic model of oxides. We show that this 
reasoning is of little validity as it assumes all atoms to be 
in the same bulk-like state. On the contrary, the addi- 
tional structural and electronic degrees of freedom at an 
oxide surface allow substantial deviations from these bulk 
properties and may thus stabilize even non-stoichiometric 
surface terminations. A similar conclusion was previ- 
ously reached also for the O-rich (0001) termination of 
corundum-structured a-Fe203 indicating that polar sur- 
faces might indeed be a more general feature of transi- 
tion metal oxides. The concentration of oxygen vacancies 
found for the polar termination of RuO2(110) at atmo- 
spheric pressures and elevated temperatures could finally 
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possibly explain the high catalytic activity reported for 
this surface. 



APPENDIX A: 



1. Vibrational contribution to the Gibbs free energy 



The vibrational contribution to the Gibbs free energy 
comprises vibrational energy and entropy, cf. eq. (15). 
Both can be calculated from the partition function of an 
iV-atomic systemEa 



Z 



3JV 

E 



dk 



E< 

n=0 



-(n+i)/3fiwi(k) 



(Al) 



where (3 = 1/kT and the u>i(k) are the 3iV vibrational 
modes. The vibrational energy is then given by 

E vih -(T,V,N) = ~lnZ, (A2) 

and the entropy is defined as 

S vih (T, V, N) = k( InZ + f3E vih - ) . (A3) 

Writing F vlh -(T, V, N) as a frequency integral including 
the phononic density of states, o-(uj), and using the rela- 
tion F vih - = E vih - - TS vih -, one arrives at 



F^-{T,u) = hu - 



kT 



1 



1 

2 + e Phu _ i 



(A4) 



^Pliuj i 



In (i _ e -l»>») 



2. Ideal gas expression for fj,o 2 (T, p) 

For an ideal gas of iV particles at constant pressure, 
p, and temperature, T, the chemical potential is simply 
given by the Gibbs free energy per atom, 



i« = 



dG 

~dN 



T,p,N 



G 

N' 



(A5) 



As the Gibbs free energy is a potential function depend- 
ing on pressure and temperature, its total derivative can 
be written as 



dG = 



dG 

df 



dT 



dG 
dp 



dp = 



(A6) 



-SdT + Vdp, 



where we have inserted the Maxwell relations for the en- 
tropy, S, and volume, V. Using the ideal gas equation 
of state, pV — NkT, the partial derivative of G(T,p) 



with respect to pressure at constant temperature is con- 
sequently 



™) =V 

dp J T 



NkT 
P 



(A7) 



In turn, a finite pressure change from p to p° results in 

G(T,p) - G(T,p°) = |* (^\ dp = (A8) 

NkTln(p/p°). 

Combining eqs. ( |A5| ) and ( |Aq) , we can finally write 
for the chemical potential of O 

fi (T,p) = l/2 f , 02 (T,p) = (A9) 

= ^o(T,p°) + 1/2 kT In (p/p°), 

which is the expression used in Section HE. 
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